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ABSTRACT 


The coagulation of dust particles under the conditions in protoplanetary disks is investi¬ 
gated. The study focuses on the repulsive electrostatic barrier against growth of charged dust 
grains. Taking into account the photoelectric effect leads to the appearance of a layer at inter¬ 
mediate heights where the dust has a close to zero charge, enabling the dust grains to grow 
efficiently. An increase in the coagulation rate comes about not only due to the lowering of 
the Coulomb barrier, but also because of the electrostatic attraction between grains of opposite 
charge due to the non-zero dispersion of the near-zero charge. Depending on the efficiency of 
mixing in the disk, the acceleration of the evolution of the dust in this layer could be important, 
both in the quasi-stationary stage of the disk evolution and during its dispersal. 

DOT 10.1134/S 1063772915070021 


1 INTRODUCTION 


The first attempts to describe the growth o f particles due to their 
coagulation were made by Smoluchowski l lSmoluchowskil[l916l) 
in 1916. He presented an equation (the so-called Smoluchowski 
equation) describing the agglomeration of colloidal particles un¬ 
dergoing Brownian motion. The methods of the coagulation the¬ 
ory were adapted to the study of processes in the late, gravity- 
assisted a ccumulation of macroscopic pre-planetary bodies by 
Safronov dSafronovll 19691) . Subsequently, these methods have also 
been widely applied to theoretical mo deling of the early stages of 
the evolution of protoplanetary disk s jSchmi tt. Henning & Muchal 
19971: l^ naka. Himeno & Id 3 l2005l: Nomura & Nakagavyal l2006l : 


Bimstiel. Dullemond & Brauerll201ol: Akimkin et al.ll2013h. 

In spite of successes in theore tical modeling of the evolu¬ 
tion of dust in protoplanetary disks iTesti et alj|2014h . many im¬ 
portant problems remain to be solved. These include so-called 
growth barriers, which hinder the formation of macroscopic dust 
particles (> 1 cm in size). The first of these barriers (arising ear¬ 
liest in the evolution of the disk) is the electrosta tic repulsion be¬ 
tween dust grains with the sam e sign of charge ( I0kuzumill2009l : 
[Matthews. Land & Hvdell2012l) . Due to the negative charge of the 
grains in deep layers of the disk, an appreciable suppression of 
the dust coagulation ca n arise, even causing thi s process to fully 
cease. As was noted in lOkuzumi et al.l ( 1201 Ibh . overcoming the 
Coulomb barrier is possible under the conditions of strong turbu¬ 
lence (a > 10^^), but this can have an adverse effect on the growth 
of larger grains due to fragmentation during collisions. Overall, 
possible ways to overcome the electrostatic barrier have been rela¬ 
tively little studied, compared to barriers to dust growth arising at 
later stages in the evolution of protoplanetary disks. 

Although the charge of dust is an important (if not the de¬ 
termining) factor in its dynamics and evolution, the coagulation 


of even neutral dust grains can be suppressed by various mech¬ 
anisms. For example, grains become more compact due to col¬ 
lisions, which increases the probability that they will bounce, 
rather than coagulate. This leads to a reduction in the growth 
of t he grains and the appearance of a so-called bouncing bar¬ 
rier dZsom et alj|201o|) . As the dust grains grow, their relative ve- 
locity increases due to their rad ial drift l Wei denschillindll977h 
and turbulence induced motion ( fVoelk et alj 19801) . If the colli¬ 
sion velocity exceeds some critical value (~ 10 m/s), the colli¬ 
sions likely lead to fragmentation, rather than ag glomeration, of 
the g r ains — the so-cal l ed fragmentation barri er dBlum & WurmI 
1200^ : IWada et ^ l2009l : iTeiser & WurmI l2009t) . Radial drift can 
remove large grains from the disk so quickly that they sim¬ 
ply are not able to grow during their mo tion toward the central 
star dBrauer, Dullemond & Hennind l2008l) . This imposes strong 
lower limits on the rate of coagulation of the particles in protoplan¬ 
etary disks at the final stages of evolution (the radial-drift barrier). 

Only collisional charging is usually consi dered in model¬ 
ing of the coagulation of charged dust grains Jo kuzumil l2009l : 
lOkuzumi et alj DPI Ibl : [Matthews. Land & Hvdd 2012h . In this 
case, under the conditions near the midplane of a protoplane¬ 
tary disk, nearly all the grains are negatively charged. However, 
in the atmosphere of the disk, where the photoelectric effect be¬ 
comes i mportant, the dust grains are instea d strongly positively 
charged IPedersen & Gomez de Castro[[201l[) . Accordingly, there 
must exist a region at intermediate heights where the grains have 
zero charge, or even opposite charges (due to the dispersion of the 
charge in the ensemble of grains), which can potentially apprecia¬ 
bly increase the coagulation rate. In our study, we investigated the 
possibility and conditions for the growth of charged dust in a layer 
of grains with near-zero charge for the typical conditions of a pro¬ 
toplanetary disk. 

Note that numerous interesting effects in dusty plasmas have 
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been discovered in experime nts under micrograv ity conditions on¬ 
board the Mir space station Iportov et alj ( l2004li . and also subse¬ 
quently on-board the International Space Station. In particular, an 
anomalous increase in the rate of coagulation of charged grains (by 
four to five order s of magnitude) was found, compared to the case 
of neutral grains dKonopka et alj|200^ . 

The structure of the paper is as follows. Section 2 describes the 
method used, and a discrete equation for the coagulation of charged 
grains with a non-zero charge dispersion is presented. Section 3 
presents the results of dust grain coagulation modeling in a proto¬ 
planetary disk for neutral dust and two cases of charged dust. Sec¬ 
tion 4 discusses dynamical aspects of the evolution of grains and 
the extent to which the adopted assumptions influence the model¬ 
ing results. The Conclusion summarizes the main conclusions of 
our study. 


2 COAGULATION THEORY APPLIED TO T HE 
GROWTH OF CHARGED GRAINS 

2.1 General Form of the Coagulation Equation 

The Smoluchowski equation can be written in two forms — in¬ 
tegral and discrete. The former considers the continuous distribu¬ 
tion of the particle masses, fijn) [cm^^ g *]. where f{m)dm is 
the number density of particles with masses within dm of m. The 
distribution f{m) is normalized to the total number density of dust 
particles = J f{m)dm. The coagulation rate for particles with 
masses mi and m 2 is specified by the coagulation kernel A(mi, m 2 ), 
which is proportional to the collisional cross section, the velocity 
of the colliding particles, and the sticking probability. In the inte¬ 
gral form, the evolution of the distribution function /(m) is given 
by the equation 

df{m,t) ^ f ft I I / i\j I 

-^- ~ 2 j ,f)A(m —m ,m jam — 

0 

-f(mj) J f(m\t)K{m,m)dm. (1) 

0 

The first term here describes the coagulation of particles with 
masses m — m' and m', leading to the formation of a particle 
with mass m. The second term is negative and corresponds to 
the disappearance of particles with mass m due to their coagula¬ 
tion with all other particles. An analytical solution to this equa¬ 
tion c an be found only for t he simplest coagulation ker nels — con¬ 
stant iSmoluchowskilligl^ and linear l lSafrono^ll969h . These ker¬ 
nels are not suitable for the conditions in a protoplanetary disk, but 
can be used to test algorithms for the solution of the Smoluchowski 
equation. A numerical approach can be used with more realistic 
kernels. A solution of this equation for th e conditions in a proto- 

S ilanet ary disk is presented, for example, in iDullemond & DominikI 

200 ^ . 

The discrete form of the coagulation equation describes the 
variations in the number densities of particles with 

masses proportional to the minimum mass of a monomer 

mo: mj. = kmQ. In this ca se, the coagulation e quation takes the 
form (see, e.g., the review ( Ichandrasekhadll943l . Chapter III, Sec¬ 
tion 6)): 

= ^ Y. NiNjKij-NkYNjKkj. ( 2 ) 

^ i+i=k ,/=l 


Note that only terms in which mi+mj = m/^ are included in the 
first (double) sum, which is equivalent to the condition i + j = k 
for the adopted mass grid. Finding an exact solution to this equa¬ 
tion requires tracing a large range of particle masses, which makes 
the grid m/^ = kmo unsuitable in practice. Therefore, a logarithmic 
grid of masses of the coagulating particles is usually used. The 
problem of conservation of mass arises, since two arbitrary masses 
from such a grid m, and mj are very unlikely to correspond to a 
k such that m^ = m, +m j. A conservative algorithm is presented 
in iKovetz & OlunJ ( 1 19691) . in which the mass m, +my distributed 
between the masses for the nearest grid points to the left and right 
mi and m^. The discrete coagulation equation can then be written 
in the form 

Nk = \YY Ni^jCukKij -NkY NjKkj, (3) 

j=\ 

i e, if njjt = max {ra„}; 

l-e, ifra,t= min {m„}; (4) 

m,;>m/+my 

0, otherwise, 

mr — (mj + mj) 

E= ^ 

nir — mi 

Even in the conservative form the coagulation equa¬ 
tion is complicated to solve. This is related to the limited length 
of the mantissa of the computer representation of a real num¬ 
ber. The coagulation of particles with masses differing by 15 or¬ 
ders of magnitude will be comparable to the rounding errors, 
even using double-precision variables in the computations. Un¬ 
der the conditions of protoplanetary disks, the range of dust 
grain masses can reach 20 orders of magnitude or more. There¬ 
fore, effective (but unwieldy) algorithms have been developed, 
that make it possible to include a wide range of masses without 
the n eed for a transition to high-precision arithmetical computa¬ 
tions iBrauer, Dullemond & HennindbOOsl) . For the goals of our 
study, it is important to preserve the simple form ([3}; therefore, we 
chose an extensive approach based on using libraries designed for 
operations with real numbers with arbitrary mantissa lengths. The 
rising computational expenses for this case were compensated by 
the application of parallel-programming methods. 

2.2 Coagulation Equation Kernel with Grain Charge 
Dispersion 

The efficiency of the coagulation process is described by the ker¬ 
nel Kjj = K{mi,mj). The case of neutral grains has been stud¬ 
ied fa ir ly well recently (see, e.g., [ Brauer, Dullemond & Hennind 
)[2()()8|);|Birnstiel, Dul lemond & Braueif ) l201olk]zsom et al.[ ) l201olk 
iBlum & WurnJ]2008l) ): 

Kij = PijH{ai + aj)^uij, ( 6 ) 

where pij is the sticking probability for grains with radii a,- and 
aj moving with relative velocity uij. We assumed that the rela¬ 
tive velocity of the grains arises only due to Brownian motion, 
which corresponds to the initial stages of dust evolution (and 
grain sizes up t o ~ 10^^ cm at distances of 1 AU; see, e.g., 
iTesti et alJ ) l2014h ). Other mechanisms determining the relative ve- 
lociti es of the grains, such as turbulence a nd radial and vertical 
drift terauer, Dullemond & Henningj 1200^ become important at 
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later stages of dust evolution than those considered here. Including 
these in models for charged dust requires a much more complex 
approach and is the topic of a separate study. 

The charges of grains in the cen t ral regions of a disk w ere 
calculated in detail in lokuzu^ (l2009ll : [okuzumi et alj j201 lal fQ). 
where the electrostatic interaction energy for large grains (with 
sizes of more than 10^^ cm) was found to exceed their kinetic en¬ 
ergy associated with collisions — the electrostatic barrier against 
dust growth referred to above. 

To analyze the influence of the grain charge on the rate of 
their growth in more detail, we considered a general for m of the 
coagu lation kernel taking into account Coulomb focusing dSnitzeil 
Il94lh : 


to io ns lead to negatively charged grains . However, as a rule (see, 
e.g.. dPedersen & Gomez de Castrol201 ih . a positive charge may be 
acquired in the atmosphere of the disk due to the efficient photo¬ 
electric effect. Therefore, there should exist a region in the disk 
at intermediate heights where the grain charges are close to zero. 
The main goal of our study was a detailed consideration of this re¬ 
gion from the point of view of the coagulation theory of charged 
particles. The main equati ons used to compute the grain charges , 
based on the treatment in iTielensI l l2005h : IPraine & SutinI (Il987h : 
IWeingartner & Drainel 1 I 2 OOI ), are presented below. 

If we consider collisions of dust grains with electrons and 
singly charged ions, the probability /(Zgr) of encountering a grain 
with charge Zgrfip can be found from the detailed balance equation 




2QiQjel 

mij{ai + aj)ujj 


Uij, 


(7) 


where Qi and Qj are the charges of the grains, = ra,my/(ra,' + 
nij) is the reduced mass, and Cp is the elementary charge. It is ob¬ 
vious that the additional condition Kij{Qi,Qj) 0 is imposed on 
Kij. This form of the coagulation kernel l[7]l is valid when the mass 
and size of the grains uniquely determine their charge. 

However, the charge of a specific grain fluctuates with time; 
in other words, there is a non-zero dispersion for the charge in an 
ensemble of grains of a given size. Let the mean charge of grain i 
be Qj and the standard devi ation of the charge about the mean be 
Oj. As a rule iTielensll2005h . the distribution of grains of a given 
size Uj over the charge q will be close to normal 




1 


o,V^ 


e 


{q-Qi?- 

2aj 


(8) 


It is not difficult to show that the coagulation kernel transforms to 
the form 


/(Zg,) [ype(Zgr) +yion(Zgr)] = /(Zg, -H)ye (Zgr + 1), (10) 

where Tpg [s^*] is the photoelectron emission rate and Jg and 
[s^*]are the accretion rates of electrons and ions. To find the total 
charge distribution function /(Zgr), the equation above is written 
for each value Zgr, which reduces to the system 


/(Zgr)=/o n 

z=Zq+1 


^pe(-^ 1) 

7e(z) 


for Zgr > Zq and 


/(Zgr)=/on' 


'7pe(z) + 7ion(t;) 


( 11 ) 


( 12 ) 


for Zgr < Zq. Here, /o = We especially note that, from a 

computational point of view, it is important to choose a value of Zq 
that is close to the maximum of /(Zgr), which can be tentatively 
estimated from the approximate equality 


^0 = / j f{p,Qj,(^j)Kij{q,p)dqdp. (9) 

The simple replacement of Kjj with K*j brings Eq.([3( for neutral 
grains to the discrete coagulation equation for charged grains with 
a non-zero charge dispersion. The sticking probability was taken 
to be unity. The process of fragmentation during collisions at high 
velocities is also important for the evolution of dust. The fragmen¬ 
tation barrier against dust growth is important for large grains; i.e., 
it becomes significant after the electrostatic barrier has been over¬ 
come. Since our aim was to search for mechanisms capable of over¬ 
coming the electrostatic barrier, we neglected grain fragmentation. 

Equation written for each k= 1,..., Amass represents a sys¬ 
tem of non-linear ordinary differential equations. We solved this 
system using an expl icit 5th-order Ru nge-Kutta method with an 
adaptive step control IPress et aljl993h . Although implicit integra¬ 
tion methods can yield a solution over a smaller number of steps, 
in practice, such methods proved slower due to the additional time 
required to compute the Jacobian. 

2.3 Computation of the Dust Charge 

The photoelectric elfect and the accretion of electrons and ions 
onto a grain are the main factors determini ng the charge of a dust 
grain under most astrophysical conditions dWeingartner & Draina 
l200lh . Thus far, only the acquisition of charge through collisions 
has been considered in relation to the dust growth in protoplane¬ 
tary disks. In this case, the high velocities of electrons compared 


■7pe(Zo)-|-7ion(Zo) ~ 7e(Zo). (13) 

Equations GD, GUl are closed with the obvious relation 

L /(Z) = l- (14) 

Z^-oo 

The accretio n rate of particles with number density n and mass 
m can be written dPraine & SutinI 1987h : 

/ 8fc7’gas 2 

7acc ~ ^ Y ^^gr7(^gr! T^gasiZg^). (15) 

Th e function / takes into acco unt Coulomb focusing. Eollow- 
ing IWeingartner & Drainel 1200ll) . we took the sticking probabili¬ 
ties for electrons and ions to be 0.5 and 1, respectively. 

The photoemission rate depends on the solid-angle- 
averaged radiation intensity Tv [erg/(cm^ s Hz ster)], the 
ab sorption cross section Cah s and the photo-ionization yield 
Y jWeingartner & Drainell200lh : 

Tpe = J Cabs(v,ri!gr)T(v,Zg]-,ag]-)^JV (16) 

An important assumption we made is that the grains are not the 
dominant charge carriers in the plasma. The depth of penetration of 
the photons and the mean free path of the electrons in a grain were 
taken to be 100 A and 10 A, respectively, and the work function of 
the photoelectrons to be 8 eV, which corresponds to silicate grains. 
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3 MODELING THE GROWTH OE CHARGED GRAINS 

To show how the charge of a grain influences its growth, the coag¬ 
ulation equation © was solved locally for each point in the disk 
for three cases. In the first, the grains were assumed to be neutral; 
i.e., their coagulation is not hindered (or facilitated) by electrostatic 
interactions. We call this case the comparison model (Model 1). In 
the second (Model 2), it is assumed that the charge of a grain is de¬ 
termined only by collisions with electrons and ions. This approx¬ 
imation is valid for dense regions in the disk near the midplane, 
where stellar radiation cannot reach. In this case, the grains be¬ 
come negatively charged, and the Coulomb interaction begins to 
hinder coagulation. Both grain-charging mechanisms are consid¬ 
ered in Model 3 — collisional and radiative. This general case is of 
the most interest. 

We emphasize that the coagulation equation was solved lo¬ 
cally and independently for each point in the disk. In other words, 
it was assumed that there was no advective term in the coagulation 
equation. This assumption is very probably not satisfied for large 
grains in protoplanetary disks, since they are subject to settling to¬ 
ward the midplane and radial drift. However, such a local approach 
is useful when searching for possible ways to enhance the growth 
of grains in some part of the disk. In this case, the distribution of 
the mean grain size over the disk visually reflects the distribution 
of coagulation rates. 

3.1 Model for the Physical Structure of the Disk 

Since we did not aim to compute the detailed physical struc¬ 
ture of the disk (density, temperature, radiation field), we chose 
the characteristic parameter s of disks around single T Tauri 
stars jWilliams & CiezallOllh . The disk was assumed to be st ation- 
ary up to an age of one to three million years (lMamaiekl2009h . with 
a surface density E, that has a power-law dependence on the dis¬ 
tance to the star I?: E = 'Lq{R/Rq)Pk We assumed hydrostatic equi¬ 
librium in the vertical direction. The radial profile of the tempera¬ 
ture was determined assuming a geometrically thi n disk and black- 
body re-radiation of the UV flux from the star (lArmitagell2O10L 
formula (2.28)). The disk was taken to be isothermal in the ver¬ 
tical direction, and the gas temperature was taken to be equal to 
the dust temperature. The stellar radiation was taken to be black- 
body radiation with the effective temperature T* and a UV excess 
with effective temperature . We chose the normalization of the 
UV excess such that The degree of 

ionization was taken to be the same throughout the disk and equal 
to Xe = 10^^, which corresponds to the ionization of carbon (we 
discuss the influence of this assumption for the model of the dust 
evolution below). The coagulation equation was solved locally for 
each point in the disk on a Vr x Nz grid. 

The influence of the assumption of stationarity and the choice 
of disk parameters on the results obtained is discussed in Section|4] 
The results of our computations for Models 1-3 with the disk pa¬ 
rameters listed in the Table 1 are presented below. 

3.2 Neutral Grains 

The case of neutral grains has been considered in many stud¬ 
ies of dust evolution in protoplanetary disks. Excluding Coulomb 
forces from consideration leads to faster grain growth for denser 
regions. Since the coagulation rate is proportional to the square of 
the particle number density, the mean size of the grains a depends 
strongly on their location in the disk (Fig. [TJ. Three million years 


Table 1. Main parameters of the model 


Parameter Value 


Ro 

1 AU 

Eo 

300 g/cm^ 

Ps 

-1 

^disk 

0.02 Mq 

t ; 

4000 K 

7’UV 

20000 K 


lo-'* 

•^mass 

64-128 

Nr 

15-30 


15-30 



10"^ 


10 “ 


10“ 


10“^ 


10“ 


10“ 


1 10 100 
R[AU] 

Figure 1. Distribution of the mean grain size over the disk for Model 1 at a 
time of three million years. 

after the start of the growth, the range of the mean grain sizes in 
the disk reaches six orders of magnitude (18 order of magnitude 
in mass). Note that fragmentation can appreciably slow the dust 
gr owth rate, thereby narro wing this range. However, as was shown 
in iwindmark et alj ( I2012) . taking into account the velocity distri¬ 
bution of the grains partially solves the problem of fragmentation. 

The detailed mass distributions of the grains f{m) for points 
A and B indicated in Fig. [T] at a radius of 10 AU and various times 
are also presented (Fig. The densities at points A and B differ 
by about three orders of magnitude. Figure [^presents the quantity 
for Model 1, which reflects the contributions of various 
grains to the total mass. The zero of time corresponds to a so-called 
MRN distribution jMathis. Rumnl & Nordsieckll977h . The growth 
of the dust at points A and B occurs similarly in terms of the pro¬ 
files of the distribution function /(m), but on different time scales. 
Figures 13 and |4] show the vertical and radial profiles of the mean 
grain size for times f = 1,10^ and 3 • 10® yr. Note that the growth of 
neutral grains is fairly uniform in the sense that Iga ~ Igf, which is 
not true for charged grains. 

3.3 Collisionally Charged Grains 

Under equilibrium conditions, the electrons have higher velocities 
than the ions, so that the grains collide more often with electrons, 
leading to a negative charge for the grains. As the grains grow, their 
negative charge also increases. Their relative velocity may also in- 
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Figure 2. Mass distribution function for the grains for Model 1 (for various 
times) at point A (upper panel) and point B (lower panel) marked in Fig.[T] 



z/R 


Figure 3. Vertical profile of the mean grain size at a distance of 10 AU from 
the star for Model 1 (for various times). 

crease. Therefore, the Coulomb factor in Eq. l[7]l can eith er increase 
or dec rease during the dust evolution. It was shown in lOkuzuinil 
( l2009h that the electrostatic interaction energy begins to exceed the 
mean kinetic energy of grains with sizes of ~ 10^^ cm or more, 
so that the growth of the grains should be appreciably suppressed. 
This threshold should depend on the ambient conditions, which can 
vary strongly over the disk. Therefore, we carried out global mod- 



R[AU] 


Figure 4. Radial profile of the mean grain size in the midplane of the disk 
for Model 1 (for various times). 



1 10 100 


R[AU] 

Figure 5. Same as Fig. [T] for Model 2. The upper boundary of the color 
scale is two orders of magnitude lower than the scale in Fig.[T] 


eling of the coagulation of (negatively) charged grains in order to 
study the conditions for the appearance of the electrostatic barrier. 

Figure 0 presents the distribution of the mean grain size over 
the disk, taking into account collisional charging. The growth of 
the dust is appreciably suppressed by its charge, and the upper 
boundary for the grain size shown by the color scale is two or¬ 
ders of magnitude lower than in Fig. [T] The main difference from 
the case of neutral dust is that the growth rate is not proportional 
to the number density. This is expressed by the fact that the grains 
grow most rapidly at intermediate radii ~ 2 — 20 AU. The grains are 
more strongly charged in inner, dense regions than in outer regions, 
making the electrostatic suppression of their growth more impor¬ 
tant there. With distance from the star, the grain charge begins to 
decrease in magnitude, speeding up the coagulation. The drop in 
density in the disk becomes important in peripheral regions, lead¬ 
ing to slowing of the coagulation. 

The mass distribution function of the charged grains also dif¬ 
fers qualitatively from the neutral case (Fig. nil. This distribution is 
bimodal; this is clearly visible for point B at a time of three mil¬ 
lion years, but is also true to a lesser extent for other times and 
locations in the disk. The origin of this bimodality is related to the 
non-geometrical relationship between the interaction cross section 
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m[g] 



Figure 6. Same as Fig.[2]for Model 2. 

and the particle size. If the mass growth rate was proportional to 
the geometrical cross section, particles with initially different sizes 
would tend to acquire more equal sizes as they grew. Additional 
factors, whether gravitational or electrostatic forces, change the ef¬ 
fective collisional cross section. This means that large particles can 
grow more rapidly than smaller ones. 

A comparison of the vertical and radial profiles of the mean 
grain size for Model 2 (Figs.|7]and[8} and Model 1 (Figs.[^and|4} 
confirms the need to take into account the dust charge when consid¬ 
ering its evolution. Note that the maximum of the profiles in Fig. [8] 
move outward during the evolution of the dust. 

3.4 General Case of Charged Dust 

The ejection of electrons from grains hy UV photons leads to ap¬ 
preciably positively charged dust in the disk atmosphere. As the 
UV flux weakens toward the midplane, the efficiency of the pho¬ 
toelectric effect decreases, and the grains become less charged. At 
some height above the midplane, the mean grain charge becomes 
zero. Upper panel of Fig. [^presents the charge vertical profile for 
grains of various sizes at a distance of 10 AU from the star. The 
mean grain charge vanishes at a height z/R ~ 0.1. Since the dis¬ 
persion of the charge is non-zero (see lower panel of Fig. lU, grains 
with charges of both signs are present in this area (point B). Thus, 
electrostatic enhancement of the grain growth will operate in the 
vicinity of lines of zero charge in the disk. 

Determining the charge of grains in the disk atmosphere is 



z/R 

Figure 7. Same as Fig.[^for Model 2. 



directly related to the radiative transfer modeling. The deeper the 
UV radiation penetrates into the disk, the deeper the layer of near¬ 
zero charge lies. This facilitates the growth of grains due to their 
higher number density. We adopted a conservative estimate for the 
radiation field, and considered only radiation from the central star, 
taking into account extinction by grains with sta ndard in t erstell ar 
parameters (in accordance with formula (5.96) of iTieleni l l2005l) ). 
A number of factors (grain growth, scattering, background radia¬ 
tion) can lead to better penetration of radiation into the disk. There¬ 
fore, the mean grain sizes in the disk atmosphere presented below 
are lower limits, and the dust coagulation rates in the disk atmo¬ 
sphere will become higher in a more detailed consideration of the 
radiative transfer. 

The difference in the densities at points A and B is approx¬ 
imately three orders of magnitude; i.e., the coagulation rates for 
the neutral case differ by six orders of magnitude. It turns out that 
this difference in coagulation rates is fully compensated due to the 
absence of Coulomb repulsion and switching to Coulomb focus¬ 
ing, and leads to an appreciable dust growth at intermediate heights 
(Fig.m. The mass distribution of the grains at point B (Fig. EJ 
is not bimodal, and resembles the distribution for neutral dust. The 
dispersion of the charge in early stages (~ 1 kyr) leads to an in¬ 
crease in the dust growth, even compared to the neutral case. The 
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Figure 9. Vertical profiles of the mean charge (upper panel) and the standard 
deviation of the charge (lower panel) at a distance of 10 AU from the star 
(for various grain sizes a). 


grains in the layer with near-zero charge are even somewhat larger 
than those in the dense midplane regions (Fig. m. 

This increase in the coagulation rate at the zero-charge line 
is a fully expected result. However, in our formulation, the mean 
charges of grains of all sizes vanish at the same zero-charge line, 
and not at various regions of space. This makes the appreciable 
thickness of the layer of active dust growth that comes about due to 
the dispersion of the dust charge an intriguing result. An additional 
contribution to the dispersion could also be made by differences in 
the work function for the emission of electrons from grains of inho¬ 
mogeneous chemical composition. The absence of an electrostatic 
barrier at intermediate heights could play a determining role for 
the overall growth of dust for a disk that is undergoing sufficiently 
strong mixing. A simultaneous description of mixing processes and 
the dust growth requires substantially more complex modeling, and 
is an interesting topic for a separate study. 


4 DISCUSSION 

4.1 Influence of Dynamical Effects 

The results presented in the previous section were obtained under 
the assumptions that the disk structure is stationary over scales of 
several million years, and that there is no drift of the dust re lative 
to the gas. Turbulence, meridional circulation iUrpinllT984h . and 
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Figure 10. Same as Fig. [T] for Model 3. The upper boundary of the color 
scale is two orders of magnitude lower than in Fig.[T] 



Figure 11. Mass distribution function for the grains at point B for Model 3 
(for various times). 
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sedimentation of the grains will lead to flows into and out of the 
layer of active dust growth. The question is whether these factors 
facilitate or hinder dust growth. The answer to this question could 
be given by appropriate modeling, but the qualitative picture can 
also be described based on general reasoning. 

On the one hand, the flow of gas through the growth layer 
reduces the residence time for a specific grain in the layer. This 
should reduce the coagulation rate, although the same material may 
pass through the growth zone multiple times in the presence of 
global vortices or circulation in the disk. On the other hand, since 
small grains are efficiently dynamically coupled to the gas, the flow 
of gas delivers small grains to the near-zero charge layer. Due to the 
lower charge inside this layer, the coupling between the dust and 
gas falls, and the grains are decelerated. They begin to grow effi¬ 
ciently in this layer, which further reduces the friction between the 
gas and the dust per unit dust mass. Thus, the residence time of the 
dust in the layer grows. In other words, a stationary flow through 
the growth layer can lead to an increase in the dust density inside 
the layer, which enhances the coagulation rate. 

As the grains grow, settling toward the midplane under the 
action of the vertical component of the gravitational force becomes 
important. The scale height, determined by the balance between 
gravitation and turbulent mixing, is smaller for larger grains. Other 
things being equal, the scale height for charged dust is larger than 
the scale height for neutral dust, due to the more effective coupling 
of charged grains with the gas. Thus, the charge of a grain increases 
the sedimentation scale height, which also enhances coagulation in 
the disk atmosphere. 

T he median lifetime of a protoplanetary disk is three million 
years dWilliams & Ciez^lTOllh . Disks around massive stars dissi¬ 
pate more rapidly than disks around solar-type stars. Therefore, the 
lifetimes for individual disks vary from one to ten million years. At 
late stages of the evolution, the dispersal of the disk occurs from 
inside out, shifting the inner edge of the disk beyond 1 AU from 
the star. In Models 1-3, the inner edge of the disk lies approxi¬ 
mately at 0.5 AU, which corresponds to the dust sublimation ra¬ 
dius. The stellar radiation is efficiently absorbed due to the high 
densities near the midplane, and does not reach distances of 1 AU, 
so that grains at height z = 0 are negatively charged. The dispersal 
of the disk means that radiation begins to penetrate into the disk 
midplane, which changes the sign of the dust charge. In this way, 
the layer of active dust growth (in the region of near-zero charge) 
moves outward along the equatorial plane, with its relatively high 
number density of dust grains. Since an increase in grain size leads 
to a decrease in opacity, this process will proceed in an accelerat¬ 
ing regime. Therefore, the growth of dust in the region of near-zero 
charge can become even more important at the stage of the disper¬ 
sal of the disk. 

4.2 Influence of Porosity and Fluctuations of the Grain 
Velocities 

The results presented here were obtained for compact, spherical 
dust grains. However, to all appearances, grains in protoplanetary 
disks have a very porous, fractal structure. This porosity increases 
the effective cross section of the grains, together with the sticking 
probability. Therefore, the presented mean grain sizes are lower 
limits for the more gen eral case of non -compact agglomerates. 
Note that it was shown in lOkuzumil (l2009h that the charge of grains 
increases with their porosity. 

The use of the mean relative velocity for grains of two spec¬ 
ified masses can distort the real pattern of the dust evolution. For 


example, considering the velocity distribution of grains of a single 
size makes it possible to lower th e fragmentation barrier against 
dust growth l lWindmark et al.l2012h . As the grain size increases, the 
relative velocity increases, and may exceed the velocity at which 
fragmentation occurs. However, grains with lower velocities will 
always be found in the ensemble, whose collisions do not lead to 
fragmentation. Consideration of the dust charge provides an addi¬ 
tional possibility for suppressing fragmentation. Since the grains in 
the midplane have the same charge sign. Coulomb repulsion slows 
particles approaching each other, lowering the fragmentation rate. 


4.3 Influence of the UV Excess and Degree of Ionization 

The magnitude of the UV excess determines how deep the radiation 
penetrates into the disk, i.e., the position of the zero-charge line. 
The more intense the UV excess, the faster the dust coagulation 
rates, since the growth layer moves deeper into the disk. 

We took the degree of ionization to be constant in the disk, 
and equal to Xg = lO^'^, corresponding to ionized carbon. The 
range of physical conditions in a protoplanetary disk is very broad, 
and th e degree of ionization varies widel y, from ~ 10^*®^ to 
unity jSemenov. Wiebe & Henning|l2004l . However, we used the 
degree of ionization only to compute the grain charges. The dust 
charge in the disk atmosphere is determined by the photoelectric 
effect, and is practically independent of Xg. The situation is the op¬ 
posite in the midplane, where collisions with electrons and ions 
determine the dust charge. However, under the conditions of elec¬ 
trical neutrality, the dependence on the number density disappears, 
and the grain charge depends only on the ratio of the masses of 
the dominant ion and the electron dTielensIlTOOSL formula (5.70)). 
Thus, the degree of ionization is important only in the transition 
layer, where carbon is ionized and Xe ~ 10^^. Note that the degree 
of ionization in deep layers of the disk could be so low that grains 
accrete all available electrons. Their charges will be smaller in mag¬ 
nitude in this case. The conclusion that the coagulation rate of the 
dust grains is substantially suppressed in the midplane of the disk is 
preserved in a self-c onsistent computation of the grain charge and 
degree of ionization dOkuzumi et al.l2011bl) . 


5 CONCLUSION 

We have considered one problem in the modern theory of planet 
formation — the electrostatic barrier against the early growth of 
dust in protoplanetary disks. Studies of the coagulation of charged 
particles have often assumed that dust grains acquire charge only 
via collisions with ions and electrons. Here, we have conducted 
modeling of grain growth taking into account radiative and colli- 
sional mechanisms for the acquisition of charge. 

Consideration of the photoelectric effect leads to the appear¬ 
ance of a line of zero charge at intermediate heights in the disk 
(where Ay ~ 1). The non-zero dispersion of the grain charge gives 
rise to a fairly thick layer in which grains with both signs of charge 
are present. This removes the electrostatic barrier in a certain layer 
of active dust growth at intermediate heights, and opens the possi¬ 
bility of removing this barrier completely for a disk with efficient 
mixing. Global flows in the disk, vortices, meridional circulation, 
and turbulence should facilitate an inflow of grains into the growth 
region from other parts of the disk. In this case, evolved dust that 
has passed through the growth region will return to other areas with 
significant electrostatic barriers. The circumstance that the scale 
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height for the sedimentation of charged grains is larger than the cor¬ 
responding scale height for neutral grains, due to the higher friction 
with gas, is important here. 

If mixing mechanisms in the disk are not sufficiently strong, 
the proposed mechanism for removing the electrostatic barrier will 
not operate efficiently in the stationary stage of evolution of the 
disk. However, the zero-charge line will pass through the disk mid¬ 
plane, where the number density of dust and the dust coagulation 
rate are high, in the stage of the dispersal of the disk. In this regime, 
the evolution of the dust can proceed in an accelerating fashion on 
short time scales, due to the positive feedback between dust growth 
and irradiation in the disk. 
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